arXiv: 1505.0105lv3 [hep-th] 12 Jun 2015 


Phase of the fermion determinant in QED3 using a gauge invariant lattice 

regularization 

Nikhil Karthilj^ and Raj amani NarayanarQ 

Department of Physics, Florida International University, Miami, FL 33199. 

(Dated: June 15, 2015) 

We use canonical formalism to study the fermion determinant in different three dimensional 
abelian gauge field backgrounds that contain non-zero magnetic and electric flux in order to under¬ 
stand the non-perturbative contributions to the parity-odd and parity-even parts of the phase. We 
show that a certain phase associated with free fermion propagation along a closed path in a mo¬ 
mentum torus is responsible for the parity anomaly in a background with non-zero electric flux. We 
consider perturbations around backgrounds with non-zero magnetic flux to understand the structure 
of the parity-breaking perturbative term at finite temperature and mass. 

PACS numbers: 11.15.-q, ll.15.Yc, 12.20.-m 


* nkarthik@fiu.edu 
1 narayanr@fiu.edu 



2 


I. INTRODUCTION 

Three dimensional Euclidean abelian gauge theory coupled to a two component massless fermion by 

Ip (A) = CTfc {dk + iAk), (1) 

can induce a parity breaking mass term for the gauge field in the form of a Chern-Simons action [in 

% tv f 

■S' = — / d^xAkFk where = tkijdiAj. (2) 

The mass term is gauge invariant in infinite volume provided the fields are assumed to vanish at infinity. This result 
has been shown in perturbation theory using the gauge invariant Pauli-Villars regularization Q. 

Theories with 2N flavors of massless fermions can have a real and positive determinant with proper pairing of 
fermions. Vacuum energy arguments show that the 0{2N) symmetry breaks down to 0{N) x 0(N) symmetry Q 
and it has been shown that dynamical masses are generated for fermions that do not break parity [g in the large N 
limit. Recent calculations along the same lines using Schwinger-Dyson equations Q attempt to identify the phase 
structure that separates one where dynamical masses are generated from others that do not. Numerical studies using 
staggered fermions [M3 have been performed and condensates have been computed for theories that do not break 
parity, again with the aim of exploring the phase structure. 

The gauge theory with one flavor of two component Dirac fermion can be regularized in a gauge invariant manner 
using the Wilson-Dirac operator [nHi3 

DUU,M) = pr.{U)-B{U)+M; ^ t (u) =-^„(U); ^^(U) = i?(U), (3) 

where U is the U(l) valued lattice link variable; .^^(U) is the naive lattice fermion operator; B{XJ) is the Wilson 
term that provides a mass of the order of cut-off for the doublers; and 0 < M < 1 is the mass in lattice units for 
the fermion. Perturbation theory computations [l3l| using Eq. m in infinite volume show that the coefficient of the 
induced Chern-Simons term in Eq. (l2])isK = — lifM>0, and k = —^ if one takes the massless limit after taking 
the continuum limit. For negative fermion masses, we would use 

D^{V,-M) = -Dl{lJ,M)=]p„{V) + B{U)-M for 0 < M < 1, (4) 

so that the induced parity breaking term as one approaches the massless limit from the positive and negative side are 
opposite in sign. 

A theory with 2N flavors with N flavors obeying Eq. ([3]) and the other N flavors obeying Eq. O can be used for a 
numerical investigation of condensates that do not break parity. We can also consider theories with non-degenerate 
fermions and arbitrary number of flavors, and study the effect of parity breaking mass terms in the limit of large 
number of flavors. 

Consider the continuum limit in a lattice simulation where we take the number of lattice points denoted by L ^ 
oo. The continuum limit needs to be taken keeping the physical spatial extent I, the fermion mass ruphys and the 
temperature T constant as L ^ oo. In a lattice calculation, it is natural to instead consider the dimensionless 
temperature t = IT and the dimensionless mass m = Imphys, measured in units of the spatial extent, to be the 
parameters of the theory and keep them constant as L —> oo. Since we study fermions on fixed gauge field backgrounds, 
the coupling constant does not play a role in the present calculations. The induced gauge action in a fixed gauge 
field background will be gauge invariant and it is of interest to study this outside perturbation theory before embarking 
on a full lattice simulation. Of particular interest is the phase of fermion determinant which contains parity violating 
terms. Consider, for example |14| . a gauge field background that has a non-zero magnetic flux. 


F^dxdy = 2 - 7073 , 


( 5 ) 


for integer ( 73 , along with a non-zero Polyakov loop, 


_ gi27r/t3 ^ associated Chern-Simons action is 


Sihs^qs) = iKTrhsQs, ( 6 ) 

and it has to remain invariant under the gauge transformation /13 —^ (/13 -|- 1). This implies k has to be an even 
integer [IM3 for this particular gauge field background in a regularization that preserves gauge invariance under 
such “large” gauge transformations [Tsl - I^ . This does not match with k = —1 01 k = —^ obtained in [13 . Effect 
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of non-vanishing gauge fields at infinity on spontaneous and anomalous breaking of parity have also been addressed 
in pTj l. In addition to the parity violating contributions to the phase of the fermion determinant, 


gir(t.m.A) ^ 

L—foo 


det Dw{U, M) 
\deiD^{U,M)Y 


( 7 ) 


in the continuum limit, there are also parity preserving contributions of the form where h(A) are integers 

associated with zero crossings of the Wilson-Dirac operator 

As a precursor to studying the three dimensional theory, consider the regularized result using Wilson-Dirac fermions 
in one dimension in comparison to the results obtained in [13, ITsl - f^ . The Wilson-Dirac fermion operator in one 
dimension is 


D^{U,M) = -1 + M + T, ( 8 ) 

where the translation operator T is {T'il>){k) = U{k)ip{k) in terms of the one dimensional link variable, U{k). The 
only physical degree of freedom is the Polyakov loop, 

L 

W = lim n U{k) = (9) 

k=l 


and the fermion determinant in the continuum, assuming L to be even, is 


( 777- \ 

u, —) 

T j / 


„i2nh 




^i2'Kh _ 


^—i2‘Kh _ 


777 = 00 

2 sin( 7 r/i)e ”^''"®5 m = 0 + 

f 2sin(7r/i)e“*'^^“*^ m = 0_ 


„—i27zh 


777 = —00. 


( 10 ) 


The result matches with the one obtained in [l9l | using zeta function regularization and has the main features discussed 
before. It is invariant under the “large” gauge transformation —>■ {h + 1) for all values of m. The part of phase 
proportional to h in the massless limit is half of its value in the infinite mass limit. As for the vacuum structure is 
concerned, the partition function for a two flavor theory with masses mi and m 2 is 


Z(mi, m 2 ) 


e mi,m 2>0 

g-(mi-m2) _|_ 2 mi > 0; m 2 < 0, 


( 11 ) 


showing that the theory with mi > 0 and m 2 < 0 is preferred over mi, m 2 > 0 . 

The aim of this paper is to study the phase r(t,m, A) in the continuum U(l) gauge field background on a three 
dimensional I x I x j torus given by 


^ 2'Kq2t 2TThi .j, . 

Ai =—^T +——+ AP; A2 = 


271(73 


2'!rh2 




A 3 = 


2'Kqit 


y -I- 2'Kh^t - 


A?, 


( 12 ) 


where qi are integers and they denote non-zero flux in the x, y and r-directions; hi S [ 0 , 1 ] denotes torons generating 
non-trivial Polyakov loops; and A^ are perturbative fields that obey periodic boundary conditions. The associated 
periodic boundary conditions on fermions are 


ip{l,y,T) = e * ! 7/((0,?/,t); 'Y{x,1,t) 


T) = e-*2’^«i^V(a;,0,T) 


-) = e~ 


‘ '0(a:,r/,O). (13) 


We refer to (73 as magnetic flux, and qi and q 2 as electric flux. The naming is not relevant if only one of the three is 
non-zero, but we also consider cases with gi ^ 0 , 92 0 and 93 = 0 in this paper and extract some results without 

completely resorting to numerical means. In addition, we numerically study the most general case with non-zero flux 
in all three directions. The phase splits into a parity even and odd part 


P = 


r 


even 


Todd; 


(14) 


with the parity even part being 


r even — 7r(qi +q2 + qs) + -K{qiq2 + q^qi + (72<73)- 


( 15 ) 
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The first term can be absorbed by changing the boundary conditions of fermions but not both the first and second 
terms. In general, the parity odd part is complicated, but it has a simple form in the case of zero temperature when 
we consider a r-dependent perturbation on a static and spatially uniform magnetic field: 

Todd = -27r/i3g3 - J (irdr'A\{ t)A2{t')G{t - r'), (16) 

where the form factor G'(r) is an odd function of t that depends on the fermion mass m and spatial torons. Our 
formulation on the lattice enables us to study G(r) without making prior assumptions concerning the local or non¬ 
local nature of the induced gauge action. We study how the form factor becomes local in the limit of m —oo and 
m —!> 0 . 

The organization of the paper is as follows. We describe lattice gauge fields on a torus in Section [TTl In Section Hill 
we derive an expression for the Wilson-Dirac fermion determinant in the lattice axial gauge allowing for non-trivial 
Polyakov loops using the canonical formalism [2^ . In Section lIVi we use the canonical formalism to study cases with 
uniform electric and magnetic fields, organized into subsections. Here, we explain the origin of the parity even phase 
in Eq. (USD. First, we present a conventional way to understand the parity breaking when there is only a non-zero 
magnetic flux. The zero crossings of the eigenvalues of the two dimensional Dirac operator are responsible for the 
parity breaking terms and the formula for the fermion determinant using lattice regularization matches the one from 
zeta function regularization We then consider the case where we have non-zero electric fluxes but zero 

magnetic flux. We show that the relevant quantity to obtain the parity even part of the phase is associated with the 
propagation of a free fermion with continuously changing momentum along a closed loop in the torus in momentum 
space in a direction defined by (< 72 ,—gi). Finally, we turn on perturbations over static magnetic field backgrounds. 
For this, we develop second order perturbation theory with in the canonical formalism in Section |V] and use it to 
study the parity odd part of the induced effective action. The results of the perturbative analysis and the numerical 
extraction of the form factor G are presented in Section IV Al and Section IV Bl 


II. GAUGE FIELD ON A TORUS 


We work on an x /3 lattice for the sake of simplicity, which can be easily generalized to a spatially anisotropic 
lattice as well. We only consider lattices where both L and (3 are even; while the continuum physics is independent of 
this choice, it helps to simplify our calculations. The spatial volume of the lattice is defined as V = L^. The spatial 
lattice points are labelled by x = {xi,X 2 ) with 1 < Xi < L, and the temporal lattice points by k with 1 < k < f3. The 
dimensionless temperature in the continuum limit is 


L 

t = lim --. 
Z/—^cx) 13 


(17) 


The continuum space and Euclidean time variables are 


and 


T = 


lim 

L—foo 


k 

L’ 


(18) 


with x,y € [0,1] and r G [O, |]. 

On this lattice, we introduce U(l) gauge fields using the gauge-links C/^(x, k). In this work, we fix the gauge such 
that the temporal gauge-links from fc = ltofc = /3 — 1 are set to identity. Non-trivial Polyakov loop variables in the 
T-direction are taken care of by the presence of C/ 3 (x, /3) = C/ 3 (x). This partial gauge fixing enables us to develop the 
canonical formalism in Section IlIIl We still have a remnant time independent gauge symmetry, g{x), under which 

Ui(x,k) ^ g\-x.)Ui{yi,k)g{yi + i) and t/ 3 (x) 5 'l'(x)C/ 3 (x) 5 (x). (19) 


In this work, we consider only the gauge fields of the form in Eq. dED. An analogous Hodge-decomposition is strictly 
true for any gauge-fields in two dimensions. In three dimensions, one should consider these as specific background 
gauge fields used in order to probe the dependence of the fermion determinant on perturbative and non-perturbative 
aspects of the gauge-field. The gauge fields in Eq. (fT^ are periodic only up to a gauge transformation with non-trivial 
winding. Since we do not require smoothness of the link variables on the lattice, such gauge fields along with fermions, 
which satisfy the boundary conditions in Eq. m, can be incorporated using gauge links and fermions that are strictly 
periodic. For our gauge choice, the lattice gauge field background corresponding to Eq. (HID is 





if xi < L 
if Xi = L, 


Ui{x,k) 
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J72(x, k) 

C^3(x) 


_ g-j^^a;i+i^^X2gi2'n-/i3 


The various background gauge fields we study in this paper are instances of the above equation. 


( 20 ) 


III. CANONICAL FORMALISM 


The partial gauge fixing defined in Section |n] naturally allows for the development of the Hamiltonian or the 
canonical formalism [2^. Let Ti{k) be the parallel transporters along spatial directions at a fixed Euclidean time, k: 


[r,-(fc)V’] (x) = 17j-(x, fc)'!/'(x +j), 


( 21 ) 


and let T 3 defined as 


[Tsijj] (x) = U 3 {K)tp{x), (22) 

be the parallel transporter that connects k = jS and fc = 1. In this gauge field background, the Wilson-Dirac fermion 
operator is 




—3 + M + — |^((Ti + 1)Ti(fc) — (di — 

^ i=l 



(d3 + l)<5'=''2-(a3-l)Tj<5'='’'5 
(da + 1 ) 5’^''’^+^ - (d3 - 1 ) (5'='’'=-! 
(d3 + l)T3J'='’l-(d3-l)5'='’^-l 


^)Thk) 


^k' ,k 


if fc = 1 


if 1 < k < P 
if k = p, 


(23) 


where the second term takes care of the periodicity in the temporal direction. In this way, we have managed to write 
the Dirac operator using operators defined on two-dimensional time-slices. The Wilson mass M is such that \M\ < 1. 
It is related to the physical mass m in the units of box length as 


m = ML. 


By using the following set of Pauli matrices, 


di = 





the Wilson-Dirac operator D can be written in the matrix form as 


D{M) = 


/-Bi 

Cl 

1 

0 

. 0 

\ 

-cl 

-Bi 

0 

0 

. 0 

Tl 

0 

0 

-B 2 

C 2 

. 0 

0 

0 

1 

-cl 

-B 2 

. 0 

0 

T 3 

0 

0 

0 

. -Bp 

Cp 

V 0 

0 

0 

0 

. -cl 

-Bg 


where 


B^=i-M-\^{T,{k)+T}{k)), 

Ck = \ (ri(fc) - Tl{k)) - 1 {r^ik) - Tl{k)). 


(24) 


(25) 


(26) 


(27) 
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Note that Bk is a positive definite operator for \M\ < 1. We closely follow in order to obtain an expression for 
the determinant of D. We first cyclically permute the columns to the left. This gives a matrix 


where 


D'{M) 


/«! 0 0 
72 a2 0 

0 73 Q!3 


0 7i>^\ 

0 0 

0 0 


0 0 0 
V 0 0 0 


ap-i 0 

7/3 oipX/ 





-BA 

-Cl)’ 





Using the formula for the determinant of the above matrix from [27|, we arrive at 



' /3 


/ 1 \ 

det D{M) = 

1 

0) 

det 

-1 

1 

_ 1 


where the hermitean transfer matrix 7k associated with propagating the fermion across the fc-th slice is 


7k — 7k — 


Bl 


-CkB-^ CkB-^Cl + Bk 


—i/^t 


The final expression for the fermion determinant is 


det 




where 


n = i- 



T. 


3 ■ 


(28) 


(29) 


(30) 


(31) 


(32) 


This is the main formula that we use repeatedly in order to understand the phase of the determinant, T, in this paper. 
Since Bj is positive definite, the phase becomes 


exp (tr) 


detI1(M) _ , „ deiTL 

I det I1(M)I ^ ^|detH|’ 


(33) 


If are the 2U eigenvalues of Y^k=p BkT^, then 


_ ^ 1 _ e. 

exp(ir) = detT 3 lilTATl- (34) 

The positivity of the hermitean transfer matrices Tk follow from the positivity of Bk since 

(A v^)Tk (^ = {u —ClvYB^^{u —civ) + v^BkV > 0. (35) 

In addition, they satisfy the unitarity property detTfe = 1. 


A. Free field theory 

In this subsection, we find the eigenvalues and eigenvectors of B (we can drop the subscript k) for free field theory, 
where all the gauge-links are set to identity. The momenta pi and p 2 in the xy-plane are 


2Tmi 


Pi = 


L 


where = 0,1 ,... L — 1. 


(36) 
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When expressed in this momentum basis, both B and C are the numbers 

2 


h= 1-M + 2^ 


i=i 


• 2 Vi ^ • • 

sm and c = i sinpi + sinp2, 


( 37 ) 


respectively. Thus T becomes 


T{ni,n2) = 


- b 


-f 


The eigenvalues of T{ni,n2) are with 


Xp = cosh 


1 + 6^ + |cp 


2b 


( 38 ) 


( 39 ) 


The corresponding normalized eigenvectors for the zero mode ( 0 , 0 ) and the doubler modes (tt, 0 ), ( 0 , tt) and (tt, tt) are 


\P+) = 


and |p—) = 


when 6 < 1; 


if & > 1 , the above \p+) and |p—) get interchanged. For other generic modes 


b±) = 


1 


\/|c|2 + (l-e±^p6)2 


±\r,h'\2 \ 1 — 


( 40 ) 


( 41 ) 


It is straightforward to extend the free theory results to a case where uniform spatial torons hi and /12 are present. 
For this, one replaces T(ni,n2) by T{ni + /ii, n2 + /i2)- 


IV. GAUGE FIELD BAGKGROUNDS WITH UNIFORM ELEGTRIC AND MAGNETIC FIELDS 

This section is devoted to gauge-field backgrounds with constant and uniform electric as well as magnetic fields; non¬ 
zero gi, 92 and 93. We first consider the case when 91 = 92 = 0 and assume that ^ 0 . This is a standard example 
to understand the role of large gauge transformation in the parity odd part of the induced action 0 m, [11 [21 
and we will show that the results using the lattice formulation are consistent with zeta function regularization. Next, 
we consider the case of static electric fields (91 0, 92 7^ 0 and 93 = 0) by reducing the problem to a free fermion 

propagation with continuously changing momentum along a closed loop in a two-dimensional momentum torus. Apart 
from providing a different perspective to the constant magnetic field case, this also leads to an understanding of a 
parity even phase 'Kqiq2- The last subsection deals with a numerical study of the general case where both the electric 
and magnetic fields are present. 


A. Uniform and static magnetic field 

Let us consider a gauge-field background with only a uniform magnetic field 93 and the toron h^. In this case, the 
matrices Tk =T are time independent. The eigenvalues of T can be written as due to its positivity. The matrix 
Ti. for this static case becomes 

and its eigenvalues 1 — are given by 

_ g±/3A)^—i27r/t3 

It is known that 93 eigenvalues of T cross unity as a function of mass when a non-zero topological charge 93 is 
present; when m > 0 , there are U -I- 93 eigenvalues , and U — 93 eigenvalues e~^ with > 0 . The determinant 
of hist expressed in terms of these eigenvalues is 

V+q3 V—qs 

detest = n Y[ 

i=i 3=1 


( 42 ) 

( 43 ) 


(44) 






FIG. 1. The continuum limit of rriQ (refer Eq. (1461) 1 at various is shown for m = 0.7 and t = 1. The spatial lattice extent 
L is specified on top of each curve. The continuum limit of obtained by 1/L extrapolation is shown by the topmost black 
line. This continuum limit matches with m = 0.7 at all h^. 


Using Eq. (EH), the phase of the determinant is 

V+q3 V—qs 

r = Trga - 27rh3g3 + ^ Imlog (l - ^ Imlog (l - ^ , ( 45 ) 

i=l j=l 


The first term irq^ is parity even, and it could be absorbed by changing /13 —/13 + i. This formula is explicitly gauge 
invariant under a large gauge transformation hs —>■ /13 + I and is a consequence of the gauge invariant regularization. At 
any hnite fermion mass, all the LXf^ have non-zero hnite continuum limits. At zero temperature, all the exponentials 
vanish leaving only the first two terms which do not depend on the fermion mass. 

In order to check consistency with the results from zeta function regularization in [isl [l^ . we computed T for 
several values oi m = ML and several values of t For the different T that we computed at finite L, the mass 

term used in the zeta function regularization would correspond to 


inciL) 

t 


= In 


tanT 

cos 27rh3 tan T — sin 27r/i3 


(46) 


We extracted the continuum limit of m,^ by fitting the results for L = 20, 22, • • • , 48, 50 using a polynomial in -h. We 
verified that the extracted value for matches with m quite well. Our checks were made in the range m € [0.2,1.0] 
and t G [0.1,3]. We show an example of the L-dependence of in Figure [U using m = 0.7 and t = 1. The continuum 
limit of rriQ is seen to match with the m used in our lattice calculation. 


B. Uniform and static electric fields 

At any finite non-zero temperature and finite volume, it is possible to have spatially uniform and static electric 
fields, i.e., non-zero qi and 52 which are integers. Also, qi and q 2 need not have the same value. In this subsection, 
we consider this case of non-zero electric fields, but with no magnetic field. The lattice gauge field background in Eq. 
EOD reduces to 


t/l(x,fc)=e*^^ 

U2(x,fc) =e-*^^ 

C/ 3 (x) = (47) 


We focus on the parity even phase arising from this configuration. At any time-slice, the above Ui and U 2 act like 
time-dependent torons hi and ^2 whose effect is to offset the momentum. Switching to momentum basis and using 
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the replacement rii ^ rii + hi, the two dimensional transfer matrix becomes 


7 ^’ =T(ni+q2^,n2-qi^ 


^ni,SiSn2,S2 ; 


( 48 ) 


where T is given by Eq. (| 38 l) for the case with torons. The n and s are the momentum indices. The product of these 
matrices is diagonal in momentum space and it is denoted as tni,n2 


k—^ 




( 49 ) 


Since T3 is already in a definite momentum (—92,71), it becomes 92 jn 2 .s 2 +<?y Thus, 


X{TkT^ 

k—0 


_ ^ ^ni—q2,sign2+qi,S2 


(50) 


We block-diagonalize the above matrix in the following way. Starting from an arbitrary momentum (ni, 712), we create 
a cycle C by moving to (ni — 92, ?t-2 + 9 i), then to (ni — 292,112 + 291) and so on till we are back at (711,77.2). This 
will occur after P steps when both P92 and Pqi become multiples of L. We refer to P as the cycle length and this is 
fixed given 91, 92 and L. The cycle C corresponds to a P x P block, and it has the following structure 


1 

l[Tk% 

_ 

k=p 

c 


/ 


t 


\t 


ni + (P-l)q 2 ,n 2 — (P-l)qi 


ni ,712 

0 

0 


0 


0 0 
711+92,771—91 0 

0 ^711+292,Til—291 


0 


0 


0 

oj 


(51) 


The full momentum space will be split into several such P x P blocks. If we choose another (771,772) that occurs in 
the above block as the initial points of the cycle, it will only permute the entries of the block and it will not change 
the determinant. Thus the full determinant of 77 factorizes into cycles with the factor from each cycle being 


[det 77 ] c = det 


p-i 

1 — ^Tii-792,712+7-91 

r =0 


which after cyclic permutation of the product of matrices becomes 

PP 


[det 77 ]^ = det 

Since products of T have unit determinant, it follows that 


1- ]qT(7ii - 92^,?T-2 + 9i ^ 
A;=0 




[det 77 ]^ = (1 — p) ( 1- 


(52) 


(53) 


(54) 


where p is the complex eigenvalue of the product of T around a cycle. The eigenvalues ^ of the full transfer matrix are 
the P-th roots of p and 1 /p in all the cycles. The complex number p characterizes the propagation along a cycle. We 
will show that there are real cycles where p is either real or a complex number with unit magnitude. If the eigenvalue 
in the real cycle switches sign as a function of mass, then it will be associated with a non-zero contribution to the 
parity even part of T. 

The product of 7 ” taken along a cycle C on the two dimensional momentum torus has the following interpretation. 
We start with some point (^, ^) on the continuum momentum torus of size 1 x 1 . We move continuously along 
the direction (—92,91) and compute the fermion propagation along a closed loop in this direction. One can formally 
convert this into an interpretation in the continuum without worrying about regularization. The integer momenta 
(771,772) cover the entire range of integers in the continuum. The continuum Hamiltonian in the r direction at a fixed 
(771,772) is 

$ (r) = (ni - 92 + 7 ) + {^2 + 9 i+ 7 ) + ctsto. 


(55) 
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FIG. 2. On the left is the real cycle TZo with coprime steps qi — = 1 on a 6^ lattice. This cycle has a length P = 6. When 

qi = q 2 ~ 3, the cycle TZq splits into 3 cycles each of length P = 2. 


We define fermion propagation as 


^(r + dr) = (56) 

in the limit of dr 0. Let ^^(oo) be the result of propagation from the vector (1,0)* and (0,1)* respectively at 
T = —oo. Then, 


[det7^](j = det [l — (^+(oo) 0 (oo))] • 


(57) 


in the continuum. 

In order to classify cycles, consider the momenta (ni, 77 , 2 ) and (L — tii, L — 772 ). From the expression for T in Eq. 

dMl), 


r 







0-3 T 


+ <72^, 





(58) 


Using this identity, we now show that if p is associated with the ( 771 , 712 ) cycle, p* is associated with the {L — ni, L — 772 ) 
cycle. Inserting cr| between the T’s in Eq. (1551) . we have 


det 


pp 

l-IlT 


k^O 


L-m -q 2 ^.L- 


"2 + ?,!) 


det 1 



(59) 


If we take the complex conjugate of the right hand side, then the product becomes a product of with the ordering 
reversed. Using the fact that T are hermitean, and after changing the variable k to jdP — k so as to recover the original 
ordering (modulo L), we obtain 


det 


pp 


1-WT [L-ui- q2^,L- 772 + 91 ^ 


fc =0 


= < det 


pp 


1 - n “ 92^,772 + 91 ^ 

k^O \ ^ ^ 


(60) 


This completes the proof. 
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We classify cycles in the following way. The cycles C and C* are conjugate if (ni,n 2 ) and {L — ni, L — 71 , 2 ) belong 
to C and C* respectively. If C and C* are the same cycle, then we call it a real cycle and denote it by TZ. If the cycle 
is real then we have 

[det = (61) 

This implies that either p is real or p is a complex number with unit magnitude. Since the eigenvalues of conjugate 
cycles are related by complex conjugation, they can be paired together to give a positive contribution to det 7^. 
Therefore, only the real cycles contribute to the phase, which can only be ±1. The following cases are possible for 
the real cycles. 

1. When p = the factor will be real and positive in the above product. 

2. When p is real and negative, the contribution will be positive. 

3. When p is real and positive, the contribution will be negative. 

Depending on the number of real cycles in the above categories, the phase could be either +1 or —1. Let us start 
by assuming that qi and (72 are coprimes. All the cycles have the same length, P, which is even for even L. Let us 

assume that (ni, n 2 ) belongs to a real cycle. Since (L — ni, L — n 2 ) belongs to the same cycle, it follows that 


(L - ni, L - 712 ) = (ni - rg 2 , ^2 +rqi) + {kiL, ^ 2 ^), (62) 

for some integers r, ki and ^ 2 - Since qi and q 2 are coprimes and L is assumed to be even, it follows that r is even and 

(ni - ^q 2 ,n 2 + = (^(1 - ki)^ (1 - ^ 2 )-^^ • (63) 

Therefore, real cycles have to contain (ni, 712 ) = (0,0), (-§,0), (O, or (-1, ■^). Let TZq denote the cycle that contains 
( 77 , 1 , 77 , 2 ) = (0,0). For every {—rq2,rqi) in this cycle there is a partner, {rq2,—rqi), in the cycle. Only (0,0), (-jjO), 
( 0 , -j) or (-j, have themselves as their partner. Since each cycle has an even number of points, we conclude that 
one of (-1,0), ( 0 , -j) or (-j, -Dalso belongs to Pq. Since the length P can have only one factor of 2, the number of 

cycles, has to be even. Since the complex cycles pair up, the two left over from (-1,0), (O, and (-j, -j) have 
to pair up and belong to another real cycle, which we call If qi and q 2 have a common factor, then we will 
assume that we choose a set of L that all have this factor while taking the continuum limit. Under such a choice, the 
common factor of qi, q 2 and L can be pulled out resulting in multiples of cycles traced using coprime steps qi and 
q2 on a smaller spatial lattice. For the sake of clarity, we demonstrate the above statement in Figure [21 for the case 

= 92 = 3 on a 6 ^ lattice. 

We now numerically show that the phase for TZq is 0 and TZji is tt for all values of qi and 52 that are coprime. This 
enables us to write the continuum formula for the phase for this case as 

Feven = 7I'(9l + 92 + 9l92), (64) 

in accordance with Eq. (I15|) . In order to maintain nnmerical stability in the computation of the product of transfer 
matrices in Eq. (ESI), we found it useful to normalize each row separately as we multiply and cumulate the normalization 
factors in a single diagonal matrix. Using this procedure we were able to work with large L and /3, thereby essentially 
seeing the behavior of cycles in the continuum limit. The top left panel in Figure |SI shows the flow of the phase from 
each cycle as a function of mass when the background electric flux is 91 = 2 and 92 = 3. The flow is close to what one 
would see in the continuum since the computations are on a 160 ^ lattice. The phase from the real cycle TZq changes 
from being tt for m < mc{L) to 0 for m > mc{L) for some positive mc{L), which becomes zero in the continuum limit 
as shown in the top right panel of Figure |31 The real cycle T^-jr has a phase of tt for all masses. The rest of the cycles 
are complex and come in pairs as is evident from the plot. The combined phase is only due to the real cycles and is 
TT for all masses above mc{L). This is consistent with Eq. (IMl) . 

It is interesting to focus on the crossing that occurs in the TZq cycle. In order to zoom in on the crossing, it is better 
to work on a coarse lattice and we picked 91 = I and 92 = 0 on a 4^ lattice and considered m G [0.188138,0.188144]. 

We look at the eigenvalue pair ^p, as the mass is changed in this very small range. The flow of eigenvalues on the 
complex plane is shown on the bottom panel of Figure [3] The eigenvalue pair starts out being positive on the low 
end of the mass region and approach unity at some rric which lies within the range. For a range of m above 777 c, P 
and i trace a jpj = 1 locus on the complex plane. Finally, the p becomes real and less than zero. The background 
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Re(p) 


FIG. 3. The flow of phase from each cycle as a function of mass for a background with fixed electric flux. Top left panel: 
The flow for all cycles with gi = 2, 52 = 3 on 160® lattice is shown as red lines. The real cycle TZo is shown by the solid blue 
line. In addition the plot shows the overall phase as a black dotted line. It is obvious through a visual inspection that the 
eigenvalues occur as complex conjugate pairs. The phase of the real cycle, and hence the overall phase, jump at an rUc (the 
points are only connected to aid the eye). Top right panel: The lattice spacing dependence of rric is shown for various values 
of qi and 52 . In the continuum limit, L —>■ 00 , the rric vanishes. Bottom panel: Flow of eigenvalues of the real cycle seen on 
the complex plane, in the region of crossing around rric. The plot corresponds to flux qi = 1, 52 = 0 on 4® lattice in the region 
0.188138 < m < 0.188144. The pair of eigenvalues ^p, flow from real negative values for m > rric to a complex pair of unit 
magnitude, and Anally to real positive values. 


with qi = 1 and the one with qs = 1 are related by a rotation. Thus, the zero crossing of an eigenvalue A'*" in the 
latter case, can now be equivalently understood as the flip in the sign of p of the cycle containing the zero mode. The 
range of m where this behavior occurs shrinks dramatically as one approaches the continuum and the value of rric 
gets closer to zero. 

As explained, when qi and q^ are not coprimes, the cycles split into N cycles each, [N = 3 in the example shown 
in Figure [2]) depending on the values of qx and 92 . Thus, all the N cycles originating from T^-o result in a phase that 
switches from tt to zero as m crosses rric from below. The other N cycles originating from always have a phase of 
TT. Thus, the total phase becomes 

r = 7r(A^mod 2). (65) 

Only when both qi and 52 are even, N can be even. Thus the expression for the phase remains as Eq. (j64[) even when 
qi and q 2 are not coprime. 

Now we proceed to add hi and /12 to the gauge field background in Eq. (SZl). The effect is to replace rii by rii + hi 
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FIG. 4. The flow of phase from all the cycles (red lines) as a function of mass on a background with fixed electric fluxes and 
torons. The left panel shows this flow for a background with flux qi = I, q 2 ~ 0 and torons hi = 0.23, /12 = 0 on 160^ lattice. 
Similarly, the right panel is for a background with flux qi = 2, q '2 = 3 and torons hi — 0.37, h 2 — 0.23. The cycle that would 
have been the real cycle TZo when torons are absent, is specially shown by the solid blue line. The overall parity odd part of 
the phase as black dashed line. The arrows indicate the expectations from Eq. (I 68 II for the infinite mass limit. 


in Eq. (15^ : 


[defHjp = det 

The full determinant still factorizes into cycles but the real cycles now become complex and the previous complex 
cycles that were complex conjugate pairs are no longer paired. If hi and /12 are multiples oiq 2 /j3 and qi/ j5 respectively, 
then it is possible to find an integer k' = k — , such that the determinant becomes 


/3P 


L - 


k=0 


Y\_T ini + hi - q 2 -,n 2 + /12 + qi- 


( 66 ) 


[detH](. = det 


pp 

i-n^' 

fc '=0 


ni 


k' 

<12-;, n2 


hiqi + ^292 
92 



(67) 


This means that at any fermion mass and temperature, the phase can only be a function of hiqi + h 2 q 2 - In the 
continuum limit, the fact that we chose a rational hi and ^2 should not matter. We proceed to compute the phase 
per cycle and the total parity odd part of the phase of the determinant numerically in order to understand the term 
in the phase that couples hi with Two sample cases studied are plotted in Figure 0] Consider the case of gi = 1 
and q 2 = 0 with hi = 0.23 and /i 2 = 0 shown on the left panel of Figure 01 This is just a rotated version of the 
case with constant magnetic flux and a temporal toron. After removing a factor of —1 from the determinant due to 
the parity even part of the phase, the parity odd part of the phase at the largest mass is consistent with —2'Khiqi as 
expected from Eq. (051) . Next, we consider the more interesting case of gi = 2, 52 = 3 with hi = 0.37 and /12 = 0.23 
shown on the right panel of Figure 01 The parity even part of the phase is again tt as in the previous case. The parity 
odd part of the phase at the largest mass is consistent with 


r = -27r(/iigi + h2q2), 


( 68 ) 


which is indicated by arrows in the plots. 


C. Uniform and static electric and magnetic fields 

Now we consider the gauge field background where electric as well as magnetic fields are present i.e., gi, g 2 and ga 
terms are all present in Eq. (gnu with no torons. We are unable to study this case analytically. Therefore, we study 
this general case by directly evaluating Eq. (EH). We check for any loss of precision by comparing the determinant of 
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91 

52 

?3 

T/tt 

0 

0 

0 

0 

0 

0 

1 

1 

0 

1 

0 

1 

0 

1 

1 

1 

1 

0 

0 

1 

1 

0 

1 

1 

1 

1 

0 

1 

1 

1 

1 

0 


TABLE I. Phase P for various combinations of q\, q 2 and 53 . The “0” represents even integers and “1” represents odd integers. 

the product of Tk to 1. Doing so, we were able to use up to 12^ lattices. We find that det D is real for any qi, and 
ga. Thus, the phase can only be ±1 and its expression must be of the form 

r = ijiTT (qi +q 2 + qs) + ?? 27 r (gig 2 + 5293 + 939 i) + r]3'^(liq2q3, (69) 

where ryi = 0 or 1. Rotational symmetry guarantees that each r]i are the same for all directions. From the last section, 
we know that = 772 = 1- We do not have any analytical argument about 773 . In Table [H we collect our observations 
about the dependence of phase on qi, 52 and 53 . The results of the table are robust and found to be the same on 
L =4, 6, 8 , 10 and 12 lattices, and for various even and odd values for the q’s. The entries with 53 = 0 reiterate our 
observations of the last subsection. The entry with even qi and 52 includes the case gi = 52 = 0, which we understand 
as due to the mismatch between the number of positive and negative eigenvalues of a two dimensional Dirac operator. 
The other cases do not offer a simple analytical explanation. The important entry is the last one where all g’s are 
odd. Since the phase is even, it implies that 773 = 0. Thus, detD has a parity even phase given by 

T = TT (gi + 52 + 53 ) + TT (gi<72 + 5253 + qsqi) ■ (70) 

V. PERTURBATION THEORY: PARITY ODD CONTRIBUTIONS 

In this section, we return to the case of uniform and static magnetic field in the presence of a uniform toron in the 
Euclidean time direction that was studied in Section ITV Al and consider perturbations and A 2 on this background. 
We expand the determinant for T-L in Eq. (j32p in powers of A^ while considering the constant flux background and 
the toron to be non-perturbative. The transfer matrix, Tk, can be expanded to second order in A^ as 

Tk = T + Tk + Sk- (71) 

The detailed expressions for Tk and Sk are given in Appendix We write 

1 

Y[Tk=T^[l + P], (72) 

k=p 


where 


/3-1 /3-1 /3-lfc-l 

P = Y, T-^-^Tk+iT'^ + Y, T-'^~^Sk+iT^ + J2J2 T-^~^Tk+iT^-^-^Ti+iT^. ( 73 ) 

k=0 k=0 k=l 1=0 

From Eq. p2L we have 


log det H — log det Hst + log det 


l+nz^PT^ 


( 74 ) 


Using standard perturbation theory in the eigenbasis of the unperturbed T, we arrive at 


log det 


1 + nz^PT^ 


= EEy 

-;_L 1 _n -L 


p±(79-l)Af-z27r^3 

fc +1 


i± k—0 


_|_ g±/3Af —i27r^3 


_ _ „±(79 — 1 )A,-—2277/13 

W e ^ “ o*±.i± 


k—0 


\ g±/3Af — 2277/13 


-"fe+l 
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3-1 


-^EEE 


„±(/3-l)Af ±(/3-l)A^-i47rh3 




j± fc=0 (l + e±/5A±-*2./«3) (l + , 


,±/3A^—227r/}-3 


^ 2 ±,j±^i±, 2 ± 


fe + 1 


EEEE—T -w-(75) 




■i±,j± -rj±,i± 


i± 


where it is implicit that the summation over i+ runs up to 1^ + qa, while that of i— up to 1^ — q^. We use this general 
second order perturbative expression to study two cases of interest. 


A. Zero temperature 


We assume that we are working away from the massless limit and therefore limi_).oo LXf are strictly greater than 
zero for all i. Since /3 = j, we see that in the limit t —>■ 0, the first three terms in Eq. GSI are real and do not 
contribute to the phase. The last term can be simplified as follows. Let 

< (76) 


-jr-i±,j± -pj±,i± 

■^k+1 ■'l+l 


where the upper-bound is independent of /3. Then, the sum is bounded above by 


^i±,j±(/3) — 


„±(/S—l)Af=—i 27 r/t 3 




1 - e 


±a±tA; 


1 _|_ g±/3A*-i2'n-/i3 ^ g±/3A*-i27r/t3 




2 (cosh (±A± ± Xf) - 1) 


j^i±,j±_ 


(77) 


Explicitly, Li+j+ and Yi-j- vanish in the /3 ^ oo limit. Therefore, we need to consider only the products of 
and terms. The phase is 

^+93^ — 93 /3 —1/c—1 

r = Trqa - 271 / 1393 -b 1 ™ E E e(^"“^(') ^ ^ Im X*+f “ (78) 

i=l 7 = 1 fc=l i=0 

The second term does not depend on /13 and therefore the contribution from the toron /13 and the perturbative part 
are independent of each other at this order. If we assume is kept finite in the infinite L limit, then we can ignore 
the second factor in the parenthesis of the second term. The term 27193/13 is independent of m and changes in multiples 
of 271 under large gauge-transformation ^ + \. Even at zero temperature, the induced gauge action is not of 

the type in Eq. (I2|) if we include fields that do not vanish at infinity. 


B. Finite temperature and hs = 0 


Our aim in this subsection is to study the purely perturbative contribution to the phase in a possibly non- 
perturbative background at finite temperature. Since we are focusing on terms of the type, we set = 0. In 

addition we only consider A^{k) and A^{k) that depend only on time. When /13 = 0, the first three terms in Eq. (ESI) 
are real even at finite /3. The phase becomes 


After writing 


/3-1 fe-l 

r = EEEE 


j-± k—1 l—O 


^±{P-k+l-l)\f±{k-l-l)\f 

^1 (l -b ) 


Im-F, 


■i±.7± 


fc-l-1 


T7i±.i± 

•'l+l 




(79) 


(80) 


we arrive at an expression for the phase, which is written concisely as 

r = -EE G(fc-/)[A?(fc)A 5 (/)-A^(/i)A?(/)]. 

k^l 1^0 


( 81 ) 
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1/L 


FIG. 5. Left panel: The form factor G(r) at fermion mass m = 10~® is shown at a temperature t = 0.25 for the free qs — 0 
background gauge-field with no torons. The different symbols correspond to different L. At all L, G(r) shows a behaviour 
at small values of r. The best fit for this power-law using the L — 3200 data are shown as the black straight line. When 
the continuum is approached by increasing L, the power-law behaviour shifts to the left on the log-log plot such that G(t) 
approaches Grir) (refer Eq. (ISTm at all r. This continuum limit Grif) is shown as the magenta curve. Right panel: The 
approach to continuum of the scaled variable [G(t) — Gr(T)] is shown at various r. At all r, it approaches 0 with a dominant 
A scaling. For very small r, there is data collapse suggesting a perfect scaling. At larger r, there are corrections to this 
scaling. However, it is the most singular t~^/L behaviour of G(r) that is important to the phase of detD. 


The form factor G is 




G{k)= e 

(i±J±) 

It satisfies the anti-symmetric property 


sinh 


(l-fc) (±A±tA±) 


4 cosh cosh 


4-Im 




Tl 


■i±J± 


(82) 


G{f3-k) = -G{k). (83) 

This expression is valid for all q 3 ,hi and /i 2 . For the free case qs = 0 that we discussed in Section Fill Al the form 
factor simplifies to 


Gik) = ^ 

P 


sinh [{P — 2k) Ap] 
2cosh2^ ” 


('j-p+.p-y j-p+.p 


(84) 


We derive the expressions for in Appendix [BJ 

The behavior of this form factor G is shown as a function of time r = in Figure [3 The observations about the 
long and short distance behaviour of G(t) seen in the two panels of the figure can be summarized in the following 
way. The G{t) has a leading j- lattice correction. However, the coefficient of shows a singular behaviour. That 
is, the approach to the continuum limit is given by 

GiTL) = Gr{T) + ^G,iT) + o(^^y (85) 


where G^ir) is the continuum limit, while the singular coefficient of the dominant j- correction is given by 


Gsir) = + O (t-^) , ( 86 ) 

for some mass and temperature dependent function /. The continuum limit, Gr{T), seems to be well described by 
the regulator independent limit obtained by replacing XpL by its p ss 0 limit, Ap i.e., 


GAt) = ^ 


m sinh [(A - 2r) A„] 
An 2cosh2[A^] 


n 


(87) 
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making use of 


A„ = lim LXp = 

L—¥0C 


+ 47r2(nf + 712 ) and 


lim Im 

L—^oc 




r; 


■p+.p- 


m 

aT’ 


( 88 ) 


for all momenta even though they only hold true for pi 0. The above observations about G are seen at all mass 
and temperature. 

Let us consider the following perturbative fields chosen such that there is a non-zero Chern-Simons action: 


^lik) = 2 sin and A^{k) = cos . 


(89) 


The phase becomes 


I-= lim ^ 


/3-1 


L—^00 lA ^^ 
fc=i 


(/3 — k)G{k) sin 


2'Kn^k 

~T~ 




( 90 ) 


where we have made a change of variable from k and I to k — I, and used the antisymmetry property of G{k) in Eq. 
(155)) . Inserting Eq. (1551) for G{k), we obtain 



lim 
L—^oo 


( 2TTn'ik\ 


.^a(fc)sin 

k=l \ P / 


lim /3 

L—>-oo 

k^l 


k^ 


■ sm 


/ 2Tm^k\ 




(91) 


The first term arising from the continuum part of G(t) can be converted to an integral. The second term that arises 
from the singular part contributes in the continuum due to the behavior. The two terms can be expressed as 

2t 

Grij) sin (27rn3rt) dr + 277713 /( 771 , t)({2). (92) 

The second term is proportional to the momentum 773 and hence it is indeed the local Chern-Simons term. It 
contributes both in the infinite mass and massless limit showing that the parity odd contribution is regulator depen¬ 
dent im. At very low but non-zero temperatures, the contribution from the first term behaves as 




7r773 


E 

,'^2—0 


m(l-e-^) 
An [A| -k 77§7r2t2] 


where rii = nt + hi, 


( 93 ) 


after integration over r. This right away makes it explicit the dependence of the phase on the torons hi and ^2 in 
the 93 = 0 background. When the torons are absent, this infinite sum suffers from an infra-red divergence when 
t ^ 0 limit is taken before the 777 0 limit. But the sum becomes zero when the two limits are interchanged. In the 

777 — >■ 00 limit, the infinite sum always vanishes. Thus, the phase from the regular term is zero in both the infinite 
and zero mass limits and only the singular part contributes to the parity odd phase in these two limits. At any finite 
and non-zero mass the contribution from the regular term is not local since it is not linear in 773. 

The above discussion shows where the parity breaking phase arises at different masses. We now present results on 
the phase directly calculated using Eq. (ED). On the left panel of Figure [51 we show the behaviour of the phase as a 
function of fermion mass, for the perturbation in Eq. (1891) on a 53 = 0 background. We show the behavior at various 
values of hi = h 2 = h, and at a temperature t = 0.1. We did the numerical calculation using lattices with L = 60, 
80, 100, 120, 140 and 160. With these, we did a continuum extrapolation for E using a fourth order polynomial in 
j-. Changing the order of the polynomial to 3 or 5 made little difference to the results. In the figure, we show these 

continuum extrapolated values. When m —oo, the phase becomes —which is consistent with a Chern-Simons 

coefficient k = —1. Using the values of phase for m < 0.1, we extrapolated the results to 777 = 0 using a fourth order 

2 

polynomial in m. These extrapolations are shown by the solid lines. The extrapolated curves smoothly approach —^ 
as 777 0 , independent of h. This corresponds to a Chern-Simons coefficient k = — which is consistent with [l 3 l| . 

At other intermediate values of m, we find a strong dependence on hi and h 2 , which is expected from the above 
discussions for the (73 = 0 case. From the right panel of Figure [U it is clear that the toron dependence of the phase 
indeed comes from Grij). As t becomes smaller, the peak gets higher and shifts to smaller values of m according to 
Eq. dMl). 
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FIG. 6. The cross-over of phase from the m —>■ 0 limit (which is F = —^) to the m —>■ oo limit (F = —^) when qs = 0. On 
the left panel, the mass dependence of F is shown for various values oi hi — h 2 = h specified by different symbols. For large 
values of mass (m > 10), the phase is — For small values of mass, the phase approaches —as seen by extrapolation (solid 
lines) using data points with m < 0.1. On the right panel, the phase F (red squares) and the phase calculated only using Grir) 
(blue line) are compared when h = 0. 



FIG. 7. The mass dependence of phase F with a flux qs = 1 in the a;j/-plane. The different symbols correspond to various 
temperatures t which ranges from 1 to |. The solid lines show the polynomial extrapolation of the phase from small m < 0.1 
to m = 0. It is seen that the phase approaches — | at large mass and the extrapolation shows that the phase approaches —j 
when m —>■ 0. 


In Figure [3 we show a similar plot for a ga = I background. We do not find any dependence on the spatial torons. 

Therefore, we show only the result with hi = h 2 = 0. The different symbols are the continuum extrapolated results 

2 2 

at four different temperatures. Using a similar procedure as in the ga = 0 case, we find the phase to be — ^ and — ^ 
in the infinite and zero mass limits respectively. At finite values of m, there is a smooth cross-over between the two 
limits. At smaller t, this cross-over occurs at smaller values of m, well described by an ^ dependence of the phase. 

The above mass dependence is clearly hi, h 2 and ga dependent. Although implicit, one could consider them as 
hi — and ga — terms in the induced gauge action, that originates from the infra-red and would not be predicted 
by a pure Chern-Simons term. 
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VI. CONCLUSIONS 

We studied the contribution to the phase of the fermion determinant in QED 3 using lattice regularization and 
Wilson fermions at finite volume and temperature. We considered non-perturbative backgrounds that contain non¬ 
zero magnetic and electric flux. In addition, our backgrounds also contained constant gauge potentials referred to 
as torons. In the absence of torons and any perturbation, we studied the parity even contribution to the phase and 
our result in Eq. dTSl) is an extension of the result in [23 - 1^ for the case of just a magnetic flux. In the presence 
of toron in the time direction and a non-zero magnetic flux, our result using lattice regularization agrees with one 
obtained using zeta function regularization [mil. We extend this result for the case with electric fluxes and torons. 
In addition to extending the result, we provide an alternate way of understanding the parity even contribution when 
one has a magnetic flux. The connection between two dimensional topology and a parity even phase is translated to a 
sign associated with the propagation of a free fermion along a closed loop in two dimensional momentum torus where 
the momentum associated with the propagation changes as one goes along the closed loop. The direction associated 
with the closed loop in the two dimensional momentum torus is proportional (—< 72 , 9 i), the fluxes associated with the 
electric field. 

The effect of finite temperature on the coefficient of the induced Chern-Simons term discussed in the past [T1 is 
addressed here. In addition we also address the issue of finite mass. We show that the contribution at zero mass 
and infinite mass only comes from the regulator but there is also a contribution from the continuum part at non-zero 
finite masses. Whereas the contribution from the regulator is local and of the Chern-Simons type with a coefficient 
that is different at zero and infinite mass [l^ . the contribution at any finite non-zero mass is not local. In addition, 
the result depends on the presence of torons in the space directions if there is no magnetic flux. This is associated 
with the eigenvalues of the free two dimensional Dirac operator depending on the torons and the eigenvalues of the 
two dimensional Dirac operator in the presence of a non-zero magnetic flux being independent of the torons [ 2 ^. 

Our studies in various non-perturbative backgrounds suggest that we can study the following class of theories using 
numerical simulation; 

N+ N- 

z= [DU] Y[[di^p[dtP+] (94) 

j=l k=l 

with 0 < and M~ < 1. The simplest one to simulate is the one that does not break parity: Set = N~ = N 

and MjT = = M. This theory with N degenerate flavors is expected to have non-zero values for fermion bilinears 

that does not break parity in the massless limit It would be interesting to perform a large N analysis on the 
lattice formalism in addition to performing a numerical simulation at small values of N. Motivated by [s^ it would be 
interesting to study the theory N~ = 0, N'^ = N and = M. In particular, one could attempt to first study this 
theory for large N semi-classically using the lattice formalism where the non-perturbative effects modify the induced 
parity odd term at finite volume and temperature away from the conventional Chern-Simons term in order to preserve 
gauge invariance. A numerical study has to address the sign problem which might be under control for large N. Since 
chiral symmetry is not relevant and gauge invariance is maintained on the lattice with Wilson fermions, numerical 
studies can be performed with the aim of studying massless fermions without the necessity to use a formalism that 
preserves chiral symmetry [H, |3l| . 
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Appendix A: Expressions for D and S 

We derive the expressions for perturbative terms J- and <S in Eq. GB). As explained in Section El we consider 
perturbative helds A^{k) which are only dependent on time k. We expand Bk and Ck to second order in perturbation 
theory 


Bk = B 


= B 




■ 


(Al) 
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Similarly for Ck'- 


i^l i^l 

= cl + cl- 

Since, only first order terms seem to contribute to the phase, we write down their expressions: 

^ ^ ^ I (^ ~i 1 


h] = 


{r^-Tf)- di = l(Ti + T/); di = -(T2 + r|). 


(A2) 

(A3) 


The Ti are the forward shift operators evaluated on a free or constant magnetic field background. Then, B ^ can be 
expanded to second order as 

3-1 


= B-^ - B-^blB-^ - B-^blB-^ + B-^blB-^blB-^. 


r = 


Tk = 


-45-' 


CB-^blB-^ 


Using the above expressions, one can trace the steps sketched in Eq. (ED to obtain 

Tfe = 7” + J'k + Sk] 

^ B-^ -B-^C'f \ 

CB-^C^ +b) ’ 

( -B-^hlB-^ -B-^cl"^+B-^blB-^C^ \ 

(bl+clB-^C^ 

-CB-^blB-^C^ + CB-^cl^^j 

( -B-^blB-^ + B-^blB-^blB-^ (- 5 - 10 ^+ 5 - 1515 - 14 ^ \ 

+B-^blB-^C^ - 5-i6i5-i6i5-iC't) 
(- 45-1 + clB-^blB-^ (bl + clB-^C^ - clB-^blB-^C^ 

+CB-^blB-^ - CB-^blB-^blB-^^ +C'5-i5i5-i6i5-iC'1' + 0 ^ 5 - 14 ^ 


V 


Sk = 


V 


-^5-145-14^ + C 5 -icf) J 


It is straight forward to obtain T and S from the above expressions in terms of bj and ci. 

Appendix B: Perturbation theory in momentnm basis 
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(A5) 


In this appendix, we derive first order terms obtained in Appendix in the momentum basis. Using the Fourier 
transforms of Eq. (IA3I) . one obtains 


Ti = 




where, 
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Using the expressions for the eigenvalues and eigenvectors of T(p), 

~p+,p_ _ ai|cp + PiC* (1 - e^pb) + P*c{l - e~^^b) + 7 i (l + - 26coshAp) 

* 4[|cp + (1 “ e^pb)"^] [|cP + (1 — e~^pb)^] 

for any generic mode. For the zero and doubler modes, it is 

\ Pi if 6 < I 


j7P+,p- _ 


I Pt if ^ > 1 - 


(Bl) 


(B2) 


(B3) 
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When p Ri 0, using Eq. (1881) . we can replace Xp with A„/L in Eq. (IB2I) for Ti and T 2 . By expanding Im(J']*J' 2 ) as a 
power series in 1/L, we obtain the expression 



(B4) 
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